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determining  the  high-altitude  electromagnetic  pulse  (EMP)  vulnerability  of  the 
commercial  Bell  Telephone  T1  Carrier  system,  and  at  developing  a  T1  Carrier 
system  specifically  engineered  to  be  EMP  hard.  The  work  described  in  this 
report  was  performed  in  support  of  these  efforts. 
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1 .  INTRODUCTION 


Investigators  have  measured  the  current  induced  by  an  electromagnetic 
pulse  (BMP),  as  produced  by  the  repetitive  BMP  simulator  (REPS),  on  a  long, 
terminated,  insulated  cable  buried  at  a  shallow  depth  near  the  air-earth 
interface.  The  measured  current  is  compared  with  analytical  predictions  of 
the  current  induced  on  the  same  cable  excited  by  an  exoatmospheric  EXP.  The 
main  objectives  of  this  paper  are  to  evaluate  the  plane-wave  electromagnetic 
field  propagation  through  a  linear,  isotropic,  and  homogeneous  conducting 
medium;  to  represent  the  BMP  coupling  to  the  buried  cable  by  the  equivalent 
distributed-source ,  lumped-parameter  network  (LPN)  model;  and  to  compare  the 
results  to  experimental  measurements  using  an  BMP  simulator. 

The  reflection  and  transmission  of  electromagnetic  waves  at  a  plane 
surface  between  the  air  and  earth  media  are  familiar  phenomena.  It  is  assumed 
in  this  study  that  the  incident  wave,  generated  in  free  space,  is  a  linearly 
polarized  plane  wave  (constant  amplitude  and  phase),  and  the  earth  boundary  is 
treated  as  a  semi-infinite,  linear,  homogeneous,  isotropic,  and  conducting 
medium. 

The  incident  electric  field  above  the  earth  is  computed  based  on  the 
measured  magnetic  field  above  ground.  Subsequently,  the  transmitted  electric 
field  below  ground  is  computed  through  the  use  of  Maxwell  equations  and 
Fresnel  coefficients.  The  procedure  used  to  calculate  the  transmitted  fields 
from  measured  magnetic  field  data  is  reported  in  a  companion  paper.1  The 
cable  under  study  is  a  1200-ft  section  of  shielded  cable  that  is  terminated  in 
a  "short  circuit"  to  ground  at  both  ends  of  the  cable.  The  cable  is  buried  18 
in.  below  the  surface  of  the  earth.  A  cross  section  of  the  shielded  cable  is 
shown  in  figure  1.  The  incident  field  is  horizontally  polarized  and  arrives 
at  the  earth  boundary  at  an  elevation  angle  of  3  degrees.  The  earth's 
parameters  of  conductivity,  o,  and  dielectric  permittivity,  e,  are  derived  as 
a  function  of  frequency  based  upon  the  universal  impedance  for  soils,  gener¬ 
ated  by  Longmire  and  Smith,2  and  an  assumed  moisture  content  of  10  percent  by 
volume. 

The  cable  has  been  arranged  in  an  arc  from  the  center  of  a  point-source 
BMP  simulator,  as  shown  in  figure  2.  Thus,  the  peak  BMP  signal  arrives  at  all 
points  along  the  cable  at  the  same  time.  The  particular  simulator  being  used 
in  the  experiment  is  the  Army's  REPS,  a  horizontally  polarized  dipole  radiator 
driven  by  a  1 -MV  repetitive  pulse  generator.  The  current  is  sensed  with  a 
Stoddard  clip-on  current  probe  (model  91550-3),  and  the  data  are  recorded  on  a 
fiber-optic  system  and  transmitted  to  an  instrument  van,  where  they  are 
converted  back  to  an  electrical  signal  and  monitored  with  a  Tektronix  7912 
oscilloscope.  The  data  are  also  digitized,  processed,  and  stored  on  disk  for 
future  use. 


1  Rolando  P.  Manriquez  and  John  F.  Sweton,  An  Indirect  Measure  of  Belov- 
Ground  Electric  Field,  Conductivity,  and  Dielectric  Constant,  Harry  Diamond 
Laboratories ,  HDL-TR-2052  (September  1984). 

2C.  L.  Longmire  and  K .  5.  Smith,  A  Universal  Impedance  for  Soils,  Mission 
Research  Corp.,  Santa  Barbara,  CA,  Contract  Mo.  DMAS001-75-C0094  (October 
1975). 
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Figure  1 .  Cross  section  of  shielded  cable 
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Figure  2.  Test  configuration  for  HOL  cable  coupling 
studies. 


The  electrical  short  circuit  at  the  ends  of  the  cable  is  achieved  by  the 
use  of  a  calcium  chloride  salt  solution  poured  into  a  3  x  3  x  3  ft  pit 
containing  a  6-ft  grounding  rod  that  is  attached  to  the  cable  sheath  by  a 
standard  telephone  press-fit  connector.  The  electric  field  incident  on  the 
cable  is  obtained  by  the  use  of  the  total  magnetic  field  measured  at  five 
points  along  the  length  of  the  cable.  An  attempt  to  monitor  the  electric 
field  at  18  in.  below  the  surface  of  the  earth  is  documented  elsewhere.1 


1  Rolando  P.  Manriquez  and  John  F.  Swat on,  An  indirect  Measure  of  Below - 
Ground  Blectric  Field,  Conductivity,  and  Dielectric  Constant,  Harry  Diamond 
Laboratories,  HDL~TR-2052  ( September  1984). 


2.  BURIED -CABLE  RESPONSE  ANALYSIS 


The  analytical  method  used  to  predict  the  response  of  the  buried  cable  is 
similar  to  the  methods  used  in  earlier  work.3'4  The  incremental  section  of 
the  equivalent  distributed-source  LPN  model  is  shown  in  figure  3.  This  model 
incorporates  frequency-dependent  passive  elements.  The  voltage,  V^,  is  the 
transmitted  E-field  for  each  incremental  section.  The  impedance  for  each 
incremental  section  consists  of  the  ground  impedance  (Zg),  cable  impedance 
(ZjJ,  and  the  inductive  reactance  ( ZjJ  of  the  insulation  gap.  The  admittance 
for  each  incremental  section  consists  of  the  capacitive  susceptance  (Yd)  of 
the  insulation  in  series  with  the  admittance  ( Yg)  of  the  ground.  Thus,  the 
transmission-line  parameters  are3 

up  juy  //26  \ 

yw)  --T  +  -5T  109  \"V/  '  (1) 

(1  +  j)T/6 

Zi(U>  "  ~  2w~T  C°th  +  j)T/6J  '  (2) 


V“>  ■  “ST  ^  (!)  '  (3) 
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where 


V"1 


2irod  u2ired 

Mlf1^)  ' 


6q  =  1 /J nfyoOg  =  skin  depth  in  the  ground, 
6C  *  1//wfM0oc  =  skin  depth  in  the  shield, 
Og  =  conductivity  of  the  ground, 
oc  *  conductivity  of  the  shield, 
ad  =  conductivity  of  the  dielectric. 


(4) 


(5) 


3E.  F.  Vance,  Coupling  to  Shielded  Cables,  John  Wile y  and  Sons,  Inc .  (1978). 
4 Michael  S .  Bushell,  Rolando  P.  Manriquez,  George  Merkel,  and  William  D . 
Scharf ,  Aurora  Test  Cell  Electron  Beam  Environment— Response  of  Large  Loop, 
IEEE  Trans.  Nucl.  Sci.  NS-30,  No.  6  (December  1983),  4558-4563. 


€g  -  dielectric  constant  of  the  ground, 

eej  *  dielectric  constant  of  the  protective  jacket  of  the  cable 
UQ  -  4w  X  10-7  H/a, 

Yq  *  1.781  ■  Euler's  constant, 
b  ■  outer  radius  of  the  insulation, 
a  =  outer  radius  of  the  cable  shield,  and 
T  *  thickness  of  the  shield. 


Figure  3.  Incremental  tt  section  of  equivalent 
distributed-source  lumped-parameter 
network  model  of  buried  cable. 


The  transmission-line  response  is  obtained  by  the  implicit  method.5  At 
each  increment  along  the  line,  the  current  is  determined  by  solving  the  matrix 
equation 

[Z]  [X]  =  [VI  ,  (6) 

where  [Z]  is  the  impedance  matrix,  and  [I]  and  [V]  are  the  vector  forms  of  the 
currents  and  voltages,  respectively.  The  impedance  matrix  is  expressed  as  a 
set  of  IZ\  coefficients  and  can  be  found  by  the  application  of  Kirchoff's 
current  and  voltage  law  to  the  equivalent  transmission  line.  The  elements  of 
the  voltage  matrix  [V]  are  the  source  terms  for  the  transmission  line,  as 
shown  in  figure  3.  Hence,  equation  (6)  is  expanded  in  the  form5 


5B.  Carnahan,  H.  A.  Luther,  and  J.  0.  Wilkes,  Applied  Numerical  Methods, 
John  Wiley  and  Sons  (1969),  440-442. 


K 


1 


I 


b1 11  +  c112 
a2X1  +  b2*2  +  C2I3 


a3I2  +  b3J3  +  C3I4 


aixi-l  +  hi1!  +  ciIi+1 


(7) 


an-1 In-2  +  bn-1 In-1  +  cn-1rn 


n-1 


anIn-1  +  bn*n 


=  V„ 


where  the  [Z]  coefficients  are 
b,  *  Z1  +  Z2  , 

C1  *  "z2  ' 
an-1  *  “z2n-2  ' 

bn-1  *  z2n-1  +  z2n-2  +  z2n  * 
cn-1  *  ”z2n  ' 
an  °  ~z2n-2  ' 

bn  =  Z2n-1  +  Z2n-2  +  Z2n  * 

Note  that  Zod<J  =  Z(ui),  Zeven  *  1/Y(w),  ^  =  E^uJ&z,  and  n  =  total  number  of 
sections. 

The  set  of  the  [Z]  coefficients  a,  b,  and  c  alone  is  called  the  tridi¬ 
agonal  matrix.  The  system  matrix  (eq  (7))  is  readily  solved  by  a  Gaussian 
elimination  method  with  a  maximum  of  three  variables  per  equation,  and  the 
solutions  can  be  expressed  very  concisely.  The  recursion  solutions  of 
equation  (7)  for  each  frequency  yield  the  currents  through  each  branch. 

The  cable  response  obtained  with  the  frequency-domain  LPN  model  used  in 
this  study  was  then  compared  to  two  solutions  based  on  (1)  Vance's  approach3 
using  constant  or  frequency-dependent  o  and  e  and  (2)  the  time-domain  LPN 


3E,  P.  Vance,  Coupling  to  Shielded  Cables,  John  Wiley  and  Sons,  Inc .  (1978). 
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approach2 * 4 *  using  constant  o  and  e.  The  number  of  incremental  sections  was 
increased  until  good  agreement  between  the  results  of  the  LPN  models  and  the 
results  of  Vance's  approach  was  achieved  for  constant  o  and  e.  The  advantages 
of  the  frequency-domain  LPN  technique  are  that  ( 1 )  a  and  e  can  be  functions  of 
frequency  and  moisture  content  and  (2)  different  E-field  values  can  be  used 
for  each  incremental  section.  Finally,  the  solution  of  the  short-circuit 
current,  Igc(w) ,  at  the  last  branch  is  inverse  Fourier  transformed6  to  yield 
the  desired  time  response,  I8C(t),  at  the  end  of  the  buried  cable. 

3.  DISCUSSION  OF  RESULTS 

The  magnetic  field,  Hx(t),  was  measured  using  an  H-field  sensor  at  a 
height  of  1  m  above  ground  and  located  on  one  of  the  faces  of  the  Stanford 
Research  Institute  cubical  sensor  box.7  The  digitized  waveforms  of  measured 
Hx(t)  at  test  points  TP1  to  TP5  (see  fig.  2)  are  shown  in  figure  4.  The 
transmitted  electric  fields,  E(t),  below  ground  at  TP1  to  TPS  are  shown  in 
figures  5  to  7  for  constant  values  of  conductivity  o  =  0.001,  0.007,  and  0.02 
mho/m  with  dielectric  constant  er  =  15.  These  figures  show  the  early-time 
variations  of  E(t)  up  to  1  ms.  Although  the  calculations  were  carried  to 
5  ps,  all  the  E(t)  beyond  1  ms  are  small.  As  shown  in  figures  5  to  7,  both 
amplitude  and  waveform  are  significantly  affected  by  changes  in  conductivity 
but  are  relatively  insensitive  to  changes  in  dielectric  permittivity. 

In  reality,  the  conductivity  and  dielectric  permittivity  are  functions  of 
frequency  and  depend  upon  the  moisture  content  of  the  soil.  Longmire  and 
Smith's  universal  formula2  is  used  to  determine  o  and  er  for  variation  of 
moisture  content.  These  are  shown  in  figures  8  and  9.  Figures  10  and  11  show 
the  effects  of  these  frequency  functions  for  10-  and  25-percent  moisture 
content,  respectively,  on  E(t)  at  TP1  and  TPS. 

The  transmitted  fields,  conductivity,  and  dielectric  permittivity  are 
input  to  the  frequency-domain  LPN  computer  program  used  to  calculate  the 
short-circuit  current  response  of  the  buried  cable.  The  impedance  and 
admittance  parameters  were  increased  by  a  factor  of  2  to  obtain  a  better 
correlation  between  the  measured  and  calculated  short-circuit  current  at  early 
times. 


2C.  L.  Longmire  and  K.  S.  Smith,  A  Universal  Impedance  for  Soils,  Mission 
Research  Corp.,  Santa  Barbara,  CA,  Contract  No,  DNAS001-75-C0094  (October 

1975). 

4 Michael  S .  Bus hell,  Rolando  P.  Manriquez,  George  Merkel,  and  William  D. 

Scharf,  Aurora  Test  Cell  Electron  Beam  Environment-Response  of  Large  Loop, 
IEEE  Trans.  Nucl.  Sci.  NS- 30 ,  No.  6  (December  1983),  4558-4563. 

6 Alfred  G.  Brandstein  and  Egon  Marx,  Numerical  Fourier  Transform,  Harry 
Diamond  Laboratories ,  HDL-TR-1748  (September  1976). 

7  B.  C .  Tupper,  R.  H.  Stehle,  and  R.  T.  Wolfram,  BMP  Instrumentation 
Development,  Stanford  Research  Institute,  report  7990,  under  contract  to  HDL, 
Contract  DAAK02-69-C-0674  (June  1972). 
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Soil  conductivity  with  varying  moisture  content 
(5  to  25  percent). 
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Figure  10.  Transmitted  electric  field  at  TP1  to  TP5  with  10-percent 
moisture  content. 


Figure  11.  Transmitted  electric  field  at  TP1  to  TPS  with  25-percent 
moisture  content. 


Figure  12  compares  the  measured  and  calculated  E(t)  with  o  and  er  varying 
with  frequency  at  5-,  1 0-,  and  25-percent  moisture  content.  Values  of  o  = 
0.007  mho/m,  er  *  15,  and  10-percent  moisture  content  yield  close  agreement 
between  the  measured  and  LPN  data  at  early  times.  These  values  were  also 
determined  empirically  from  correlation  between  the  measured  and  calculated 
transmitted  electric  fields  below  ground,  as  cited  elsewhere.1  The  comparison 
between  the  measured  and  calculated  Isc(t)  with  o  *  0.001,  0.007,  and  0.02 
mho/m  and  er  *  15  are  shown  in  figure  13.  Figure  14  shows  the  late-time 
measured  short-circuit  current. 
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Figure  12.  Comparison  between  measured 
and  LPN  data  with  varying 
moisture  content  (o  and  ef 
both  functions  of  frequency). 
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1  Rolando  P.  Manriquez  and  John  F.  Sweton,  An  Indirect  Measure  of  Below- 
Ground  Electric  Field ,  Conductivity,  and  Dielectric  Constant,  Harry  Diamond 
Laboratories,  HDL-TR-2052  (September  1984). 
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Figure  14.  Late-time  measured  short-circuit 
current. 


4.  CONCLUSION 

It  can  be  seen  from  the  data  (fig.  12)  that  the  risetime  and  amplitude  of 
the  calculated  and  measured  currents  are  in  good  agreement.  However,  for  the 
late-time  waveshape  (fig.  12  and  14),  there  is  a  significant  divergence 
between  the  measured  and  calculated  data.  The  reason  or  reasons  for  this 
disagreement  are  not  known,  but  it  is  clear  that  if  this  discrepancy  results 
from  the  (near-field)  proximity  of  the  cable  to  the  simulator  (a  likely 
possibility),  then  the  result  to  distributed  system  testing  could  be 
significant.  On  the  other  hand,  if  this  discrepancy  is  the  result  of  some 
factor  not  adequately  accounted  for  in  the  analytic  calculation,  the  result  to 
distributed  system  analysis  would  be  similarly  significant. 
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